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Abstract 

Observed optical reflectivity in the infrared spectral region is compared with 
theoretical predictions in a strongly coupled electron-phonon system. Start- 
ing from a Frohlich Hamiltonian, the spectral functions and their temperature 
dependence are derived. A full analysis including vertex corrections leads to 
an expression for the optical conductivity cr{uj) which can be formulated in 
terms of the well known optical conductivity for a quasi-isotropic system with- 
out vertex corrections. A numerical comparison between the full result and 
the so-called "extended" Drude formula, its weak coupling expansion, show 
little difference over a wide range of coupling constants. Normal state optical 
spectra for the high-Tc superconductors YBa2Cu307 and La2~xSra;Cu04 at 
optimal doping are compared with the results of model calculations. Taking 
the plasma frequency and eoo from band structure calculations, the model has 
only one free parameter, the electron-phonon coupling constant A. In both 
materials the overall behaviour of the reflectivity can be well accounted for 
over a wide frequency range. Systematic differences exist only in the mid- 
infrared region. They become more pronounced with increasing frequency, 
which indicates that a detailed model for the optical response should include 
temperature dependent mid-infrared bands. 
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I. INTRODUCTION 



The optical spectra of metals in the infrared (IR) spectral region depend sensitively on 
the interaction between electrons and phonons. Deviations from the theoretical spectrum 
without any phonon contribution stem from the so-called Holstein mechanism in which 
the incident photon is absorbed in a second-order process involving creation of both a phonon 
and an electron-hole pair. The detailed description of this mechanism was given by Allen 
. Despite the fact that the physical mechanism of electron-phonon coupling is rather well 
understood for almost three decades, little process was made in the systematic experimental 
study of coupling effects in optical (IR) spectra 0, nor has the phenomenology been explored 
beyond the lowest order effects. The lack of systematic experimental investigations appears 
to be related to the following two reasons: 

Firstly, historically measurements of optical spectra of metals where limited to photon 
energies uj ^ 0.05 eV ||^. We shall argue in this paper that the effects of the electron-phonon 
interaction on the optical spectra of ordinary metals are rather weak at these energies. The 
optical conductivity in this spectral range can usually be described in the framework of the 
Drude formula 

a{u) = ^ (1) 

where 1/r is the relaxation rate of electrons due to their interaction with impurities and 
phonons. It can be calculated using the commonly adapted Bloch-Griineisen-type formula. 

Secondly, the measurements of the optical conductivity are complicated in ordinary met- 
als by the anomalous skin effect which is mostly present at low temperatures and in the 
far-infrared (FIR) spectral region. It leads to serious difficulties in the interpretation of 
reflectivity and absorbtivity measurements. Only a few observations [^,0 are known to us 
where Holstein processes where identified in the normal state of ordinary metals. 

The discovery of the high-T^ superconductors (HTSC) has dramatically changed the 
experimental situation. First of all, the experimental methods were improved radically by 
extending the accessible energy range down to ~ 10 cm~^ and by increasing the accuracy 
of the measurements. Secondly, it appears that HTSC systems allow the observation of the 
electron-phonon interaction in the optical spectra more easily and more clearly than it was 
possible in the ordinary metals. Thirdly, there is no anomalous skin effect in these systems 
for light with electric field parallel to the Cu-0 planes. 

It is the purpose of this paper to review the theoretical situation. We then extend the 
treatment of the Frohlich Hamiltonian to include vortex corrections. Theoretical predictions 
are then compared with experimental observations following earlier attempts to connect 
the IR reflectivity and absorption spectra of HTSCs with features of strongly interacting 
electron-phonon systems 

We shall show that most but not all observations can be described rather well within 
such a scheme and identify some pertinent open questions. 
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II. DERIVATION OF THE OPTICAL CONDUCTIVITY IN THE FRAMEWORK 

OF THE FROHLICH HAMILTONIAN 



We start from the description of metals with electron-phonon interaction by the standard 
Frohhch Hamiltonian 

H = J2 ^Ki Ck,i + J2 ^k(q, h i', A) cl i Ck+q,i' X (2) 

k,i k,q,i,j',A 
q,A 



X 



Here the first term is the kinetic energy of an electron with given momentum k and band 
index i, the last term is the energy of the phonon with momentum q and mode A. The 
second term represents the electron-phonon interaction, where gk{c[, i, i', A) is the matrix 
element of this interaction. The use of this Hamiltonian for the self-consistent calculation 
of the electron and phonon Green's functions cannot be rigourously justified in the general 
case. It was shown p|-p!l| that this Hamiltonian can be used for the calculation of the 
influence of the electron-phonon interaction on the electronic properties for systems where 
no low-energy collective excitations of electrons are present. First-principle calculations 
10| of the physical properties of a number of ordinary metals have shown that the Frohlich 



Hamiltonian is a very good starting point for the analysis of all features related to electron- 
phonon interactions. In HTSC materials we can not expect this Hamiltonian to describe 
the optical response completely because additional low energy excitations (e.g. spin and/or 
charge excitations) occur. However, these excitations contribute little to the general trend 
of the "anomalous" optical properties of HTSC systems at high energies and are simply 
ignored in this study. The question we wish to answer here is the following: to what extend 
can observed optical spectra in the normal state be understood in terms of the most simple 
model of electron-phonon interaction as represented by equation (^. 

The many-body electron-phonon interaction are usually calculated by Green's function 
method |]T2[|. Let us introduce the electron and phonon one-particle thermodynamic Green's 
functions 

G,(k,r) = -(T.4,(r)ck,(0)) (3) 

and 

Z},(q,r) = -(T.6i,(r)6q,(0)). (4) 

The Wick operator T,- reorders the operators following it in such a way that r increases from 
right to left. For non-interacting particles the Fourier components of the Green's functions, 
Gi{k,iun) and Dx{q,iuu), have the very simple form 

G°(k,i^„) = : (5) 



and 



^^(q,i^.)= (-^ — ] ■ (6) 



3 



Here iujn = (2n + l)7rT and iuji, = 2z/7rT are the Matsubara frequencies for fermions and 
bosons, respectively, and T is the temperature. The value ek,j is the electron band energy 
and, correspondingly, u^x is the phonon energy for the Ath mode. In the following we do 
not consider the renormalisation of the phonon Green's function due to the electron-phonon 
interaction. For convenience we present the phonon function D\{q, icjj,) in the spectral form 

oo 

D,{q, icu,) = - fdnim D,{q, Q) (-^-^ - -^Tt) ■ (7) 

TT J \lUJu — i'i 1UJ + \1/ 



For the non-interacting case the spectral density lmDx{q,Q) has the form 

lmDl{q,n)=7T6{n-uj^x) ■ (8) 

In the following we tread the spectral function lmDx{q,^) as an experimental quantity 
and calculate the influence of the electron-phonon interaction on the electronic properties. 
These properties are reflected in the one-particle electron Green's function and the optical 
conductivity aa/siuj). The one-particle Green's function for electrons in the presence of the 
electron-phonon interaction |]13| can be written as 



G-\k,iu;n) = G'o'(k,i^„) - S(k,ia;„) (9) 

where S(k, ia;„) is the electron self-energy part. One of the main results of Migdal |]13| is 
that the electron self-energy can be calculated using the simplest first-order approximation 
in the electron-phonon interaction, neglecting all vertex corrections as being small of the 
order of cud/cf- Here ub is a characteristic phonon energy and ep is the Fermi energy of the 
electrons. An analytical expression of S(k, icUn) is 

Ei{k,iuJn) = -Tj2 E |^7k(k-k',z,z',A)|'Z}A(k-k',ia;,) x (10) 

(^v k',j',A 

xG(k', iujn - iuj^) . 
The summation over the momentum k' is represented in integral form 

oo 

E = E/ deE^(e-^M) (11) 

'^1* * — oo ^ 

so that Sj(k, iu;„) becomes 

oo 

E,(k,icu„) = -T^ f de |^k(k-k',z,«',A)|'5(e-ek',.') X (12) 

oo 

X Jdnim Dx{k - k', n) G(k', iu;„ - iu^) x 



X 



iuj,, — Q luo,, + VL 



The analysis of (p^, as performed by Migdal, shows that the essential values of iun and iu^ 
are of order wd- This means that small frequencies of order uj-q are also significant for e. In 
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this case e can be neglected in S{e — ek,j). The electron Green's function is used in the form 
given by Eq. (||) leading to 



oo 

s,(k,icu„) = -Tj2Y.^M I a,'(k,k^^])F(^]) i - 



1^1/ k' 
oo 

X / de 



(13) 



where the spectral function of the electron-phonon interaction 



a: 



■(k, k', n)F{n) = |(7k(k - k', I, i\ A) I' Im DA(k - k', n) 



i'X 



(14) 



was introduced. The sum over the cJj, = 27rz/T can be easily performed with the result 

1 1 \ 1 



(15) 



_ N{n) + 1 - fje) ^ Njn) + f{e) 
iuin — r2 — e iujn + — e 

where N{Q) and /(e) are the Bose and Fermi function, respectively. To obtain the one- 
particle Green's function describing the electron excitation spectrum we use the analytic 
continuation of the self-energy on the imaginary axis u. This can easily be done by changing 
iuJn in Eq. (^) to ut. Consequently, the final expression for the self-energy reads 



oo 



(16) 



where 



L{uj ± i6, n)= J de 



Ar(^]) + 1 - /(e) ^ iV(^]) + /(6) 



u — Q — e ±i6 Lu + Q — e ±i6 



(17) 



The function E^('^)(k, to) denotes the retarded (advanced) self-energy which is an analytical 
function of the variable u in the upper (lower) half of the complex plane. The integral in 
(p!7|) can be evaluated analytically, yielding 



L{uj,n) = -2ni 



N{n) 



UJ 



2nT 



1 .Q + uj 



27rT 



(18) 



where "^{z) is the digamma function. The self-energy Sf'^(k, w) expressed by Eq. (|T6D 
depends only on the direction of the momentum k on the Fermi surface. It is convenient 
to present this dependence by expanding all functions involved in this expression over the 



complete and orthonormal set of functions introduced by Allen |T^. These so-called "Fermi 
surface harmonics", ^E'j(k), satisfy the condition 
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Y: ^,(k)vl/,v(k) 5(6k - 6) = S,,,N{e) , (19) 

k 

where 

Nie) = J26{e^-e). (20) 

k 

In terms of this set we write 

Sf'^(k,.;) = Es5^(^)M/,(k) (21) 

j 

Gf'\k,u;) = j:Gf:f{e^,u;)^,{k) (22) 

j 

a^,(k, n)F{n) = E E {la^'i^ - k', h A)|' lmDx{k - k', ^)} . x (23) 
xv[/^(k)v[/,v(k') . 

It follows from Eqs. (|2l| - ^3D that the self-energy coefficients E^j^(a;) are given by 

Effiiu) = E / alj,{n)F{n)L{LU ± 15, r^) , (24) 

y 

where 

4,,(n)F(n) = ^ X:E^(ek)5(ekO«^(k,k',f^)F((^)vI/,(k)^,.(k') . (25) 

Here A^(0) denotes the density of electron states on the Fermi surface. The first Fermi 
harmonic is ^E'olk) = 1 and the function 

alo{Q)F{Q) = ^ E 5(6k)5(6kO«^(k, k', Q)F{Q) (26) 

is obviously the well known Eliashberg spectral function which determines the superconduc- 
tivity of metals in the simple s-pairing case. Introducing the real and imaginary parts of 
the self-energy, Si(k, cj) and S2(k, cj), respectively, the Green's function becomes 

^-^(k, + 15) = - ek - Si(k, + i6) - iS2(k, uj + i6) . (27) 

The pole of the Green's function determines the spectrum of one-particle excitations. At 
small energies Eq. ( P?]) can be rewritten as 



9Si(k,c<;) 



du 



iS2(k,cj + i5) . (28) 



Then the pole of G occurs at ujq, which is given by 
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ujq = E] 



1 

2^ 




S2(k,Ek) . 



Here 



Aw = - 



aEi(k,o;) 



duo 



(29) 
(30) 
(31) 

(32) 



1^=0 



describes the change of the effective mass of the electron, while ^ describes their relaxation 
rate. All functions can be rewritten in terms of Fermi harmonics and the spectral function 
a|j/(f2)F(n). We shall turn to this problem later when we discuss the conductivity of metals 
in the presence of electron-phonon interaction. 

As usual we write the conductivity of a metal in the presence of electron-phonon inter- 
action in terms of the analytically continued electromagnetic kernel Kap{oj), 



Attiuj 



(33) 



which, in turn, is expressed through the one-particle Green's function Gi(k,uj) and cor- 
responding vertex function Tp. In the framework of the thermodynamical theory of the 
perturbations the expression for Ka/3{uj) is 



K. 



0/3 



T E E icu„OG'(k', icOn' + ico) x 

xr^(k',iu;^/,iu;„/ + iu) , 



(34) 



where v^, is the a-component of the electron velocity. Using the ladder diagram approxi- 
mation the equation for the vertex function is written as |^ 



T/3{'k',iUn',iuJn' + iujn) = v^, + T E Kfl'k"(k' - k",i,^', A))|^ X 



(35) 



k",n" 



xDx(k' — k", iun' — iuJn")G(k", ico'„//)G(k", iun" + ii^n) x 
xTaCk", iun", iu;„" + iu;„) . 



Before we solve Eqs. ( p^D and (p5D we simplify them somewhat. Firstly, as the conductivity 
of any metal in the absence of a magnetic field can be diagonalised in the appropriate 
representation we omit in the following the indices a and f3 considering the conductivity as a 
scalar. This is done bearing in mind that the absolute value of the conductivity in non-cubic 
crystals is anisotropic and that the corresponding functions determining the temperature 
and frequency dependence of a{u,T) reflect this anisotropy. Secondly, we omit the electron 
band indexes taking into account that interband transitions can be calculated separately if 
needed. We also omit the electron spin index multiplying by two the sum over the electron 
momentum k. After that Eq. (p5|) is rewritten in a simplified form 
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oo 

r^(k', iujn', Wn' + \ujn) = t;^/ + 2T ^ ^ vl,„ / de 8{t - Ck") x (36) 

oo 

X- / dQa\k',k")F{Q) X 

TT J 


^ 1 1 A 

y f I y 

xG(k", iti;„")G(k", ic<j„// + icUn) x 
xrj;(k", iu;„", ico-n" + iujn) ■ 

To establish some important steps in the derivation of the general formula for (j{uj) we 
consider, as a first step, the case where the vertex correction to the bare vertex can be 
neglected. Then the expression for the electromagnetic response kernel K{iujn) becomes 

oo 

K,^{iu,^) = 2i:Y.Yl I de5(e-ek")(^k')'x (37) 
xG(k', ict;„/)G(k', iun' + loJn) ■ 
We can now use the Poisson summation formula 

E n-^n) = I do; ^ ' (38) 

n=-oo 27riT J eT + 1 

where the contour C encircles the imaginary cj-axis. After that we expand the tu-contour 
to infinity, picking up contributions from the singularities of our integrand at ia;„ = ek' 
and iuJn = Ck' — ii^n- As a result we find after rather lengthy but simple calculation the 
analytically continued electromagnetic kernel as 



oo 



K{uj) = 2Y,{vl,f I de5(e-ek") / du' x (39) 

k' 

— OO — oo 



X I tanii I — I — tanii I — — ^ | | x 

xnj;^^(k',c^',cu) . 

Here n^^(k', cj', cj) has the form 

n^^(k', cu', uo) = G^(k', uo' + uj)G\k', uo') , (40) 

where ^^(k',^^' + oj) and G^ik'^uj') are the retarded and advanced Green's function, re- 
spectively. Their self-energy parts are given by Eqs. (|1^ - [1^). Just as in the case of the 
one-particle Green's function the relevant values of uj are small in comparison to the Fermi 
energy, the value of e in 5(e — ek) can be neglected, and the integration over e can be carried 
out. As result we find for the conductivity 
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aiuj 



n<m..)_l cl.' (tanh (^) - tanh [^)) x (41) 

where the function no(kF, uj', uj) depends on the position of the momentum kp on the Fermi 
surface and is given by 

u + i;^(kF, UJ + uj'j - S^(kF, uj') 

For the quasi-isotropic case we get, of course, the well known result for the optical 
conductivity 



aiuj) = ^ J duj' (^tanh (^^^) - tanh (^^^ ) no(c., u;') , (43) 

where the plasma frequency of electrons ujp\ is 

^?i = 2e^EK)'5(ek) (44) 



UJ„ 



and the function no(ti;,ti;') is 

Here we introduced the relaxation rate from impurity scattering ^— . 

'Hmp 

For the anisotropic case we should expand all functions under the integral in Eq. (^if ) 
over the Fermi harmonics. The value (v^)"^ is just the square of the Fermi harmonic of the 
order = 1 (for details see Ref. 0]). The expansion of the function IlQ{kp,u',u!) can be 
written as 

no(kF,^',c^) = J2^,iuj',uj)F^ikF) . (46) 

j 

The non-zero result for the conductivity will arise from the first harmonic, FjCkp) = 1, and 
from the harmonics with even order N > 2. The first harmonic gives the same result as 
found for the isotropic case. An example of a higher harmonic is the one which transforms 
as the representation of the crystal symmetry. It has the form 

= ''y--'''y (47) 

//o o\o\ ' 



At this point we are not certain to which extent the higher harmonics are relevant in HTSC 
systems but it is known that their influence is small in normal metals. This question should 
be considered in more detail but is beyond the scope of this paper. 
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While the analytic continuation of Eqs. (^), ( p4D in the zero-order approximation and 
absence of the vertex function r(k, icu^, icj^) is straightforward, it becomes a non-trivial 
task in the presence of r(k, icj^, itUm)- The difficulty arises from the existence of a manifold 
of functions which can be obtained as a result of analytical continuation for one variable 
while another is fixed. However, this difficulty can be avoided by changing, as above, the 
sum over the Matsubara frequencies to a contour integral. The contour consists of (three) 
circuits around the imaginary axis of the variable cj', avoiding all poles and branch cuts of 
the integrand. Using the methods developed in Refs. [|l],|16| we obtain 



oo oo , , , 



2T 



tanh 



UJ 

2T 



X 



(48) 



xn^^(k,o;',o;)r,(k,o;',o;) 



and 



oo oo 

r,{k,u',u)=v^ + 2Y^ f de(5(e-ek') / 



du" 

27lT 



(49) 



X 



tanh 



2T 



UJ" 



Akk'(t^' - + W) - tanh — Akk'(t^" - uj' - i6) + 



2T, 



io'-io" 



+coth , _ 
V 2T 



(Akk'(^" - ^' + i^) - Akk'(c^' - oo" - \S)) 



where the function 



oo 

Akk'(^) = / dfi a^(k,k',fi)F(fi) 



1 



1 



— + f2 



(50) 



was introduced. Taking into account that apart from no(k', a;", a;) all functions under the 
integral on the right-hand side of Eq. (^) depend only weakly on the variable e\j we can 
integrate over this variable and find 



r,(kF,a;',a;) =i;^ + ^5(ek0 / d^" no(k^, ^) x (51) 

k' 

— OO 

X \i{uj' - i(5, uj") - i{uj' + uj + is, n, oo")] r^.(k^, u", oo) . 

Here Ho(kp, cu", cj) is defined by Eq. (]4^ ) and I{oo,Q,oo') is 

^''^ cj-fi-cu' oo + n-oo' ^ ^ 

At this stage it is useful to introduce a new function 'jxi^.F , oo' , oo) , 

^x(Mf, u', u) = Ho(kF, uj', u)Tx(kF, u', u) . (53) 
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In terms of this function the conductivity can be expressed as 



X7a.(kF,u;', . 
For 7a;(kF, cc'', one can write the equation 

oo 

cu7^.(kF,cu',w) =vl + Y.5{e^>) I dfia2(k,k',fi)F(fi) x (55) 

k' 

oo 

X j <\u"[I{uj' -i5,Vt,uj")-I{u' + uj + i5,n,uj")]x 

X (7x(kF, u) - 7x(kF, ^\ ^)) ■ 
Now we use the expansion of the functions in Eq. (|55D in Fermi harmonics. This yields 

oo 

u^^iuj', uj) = Ivfj"'^ 5,, + fdQ X (56) 

r 

oo 

X J duj"[I{uj' -i6,n,uj")-I{uj' + uj + iS,n,uj")]x 

— oo 

X ia].,{n)F{n)-fjiuj",uj) -Y.Cn'j"a%{n)F{n)-fjiuj",uj) 

\ j" ) 

where we used the Clebsh-Gordon coefficients 

= ]^ E^(ek)vl/,(k)vl/,v(k)vl/,„(k) . (57) 

All following calculations are greatly simplified if 

4,(fi)F(fi) = ^X:E^(ek)5(ek')«'(k,k',l])F((])v^,(k)v^,,(k') (58) 

has the diagonal form 

a]^,{^)F{^) = a]in)F{n)6jj, . (59) 

In ordinary metals this assumption is well satisfied, where a^(k, k', Q)F{Q) depends mainly 
on the difference of the momenta k— k'. There are also some restrictions on the non-diagonal 
matrix elements a'jj,{Q)F{Q), j ^ j', connected with the crystal symmetries (see e.g. Ref. 
1^). Nevertheless, it is difficult to say something definite about the accuracy of this as- 



sumption in HTSC systems without concrete calculations of the functions a^(k, k', fl)F(fl). 
Such calculations where never made until now. We use this approximation as a first step. 
In addition, we restrict ourself to the use of the first two Fermi harmonics. 
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It is useful to search for the solution of Eq. (|56D in the form 

1 



(60) 



where, after lengthy but simple calculations, the equations for the transport "self-energies" 
can be written in the form 



oo oo 



SS(a;' + a;,a;') = ^EE /^^ / d^" + a; + 15, fi, l^7k(k' - k, A)|' x (61) 



k k',A 



-oo 



xImZ}A(k-k',u;')^ I 1 



and 



oo oo 



^) = ^ E E / dfi / du" I{co' - iS, n, CO") \glik' - k. A) 



X 



(62) 



xImDA(k-k',cu' 



1 - 



v^uj + St^fcu" + 00, io') - J:t(io'' iO') 



The Eqs. (|6l|) and (|6^) are still implicitly integral equations for E^J'^. However, the integrand 
on the right hand side of Eqs. (0)and ( |62D depends only weakly on the functions S^'^. The 
detailed numerical analysis of these equations, given by Allen for the case T = 0, has 
shown that the assumption 



(63) 



lj-J:1{lo' + uj,lo')-EUuj',^) ^ 
CO - Sg(cj" + CO, CO') - S4(cj", to') 

is satisfied with good accuracy. The only difference in Eq. (|6TD compared to Eqs. (^Sf ) 
and (^Sf ) is the appearance of the transport "self-energies" instead of the one-particle 
self-energies. Accordingly, the equations for E^'^ can be written in a form which largely 
resembles the equations for the one-particle self-energies, namely 

oo oo 

El{oo' + oo,io) = Jdn al{n)F{n) J dio" I{io' + io + i6, fi, oo") , 

-oo 
oo oo 

JdQal{Q)F{Q) J du" I{oo' -iS,Q,uj") , 



^tAu\uo) 



(64) 
(65) 



and 



oo 

al{^)F{^) = j^^Y.Y.5{e^)5{e^.) j duo" |<7^(k - k', A)f ImZ},(^]) x 



k k',A 



X 



1-^1 . 



(66) 
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The expression for the conductivity can then be written in a form which is formally identical 
to Eq. (^31) obtained without vertex corrections 



-oo 



00 + EfJuj' + UJ)- Etriuj' + UJ) + ^ ■ 

"Hmp 

Equation ( |67D has just the form which is normally used for the calculation of the transport 
properties of metals (see. e.g. |llO| , [T^ ). 

The expression ( |67D for the conductivity can be simplified even more in the case of weak 
electron-phonon interaction 0, to yield the so-called "extended" Drude formula 

^(^) = ^ .J. , (68) 



in iu - W{u) - ^ 



where 



= iJi 1 , (69) 

\^ m J Ttr{uj) 

oo 

Wiiv) = -2i I d^«u-(fi)m^(^,^) , (70) 



and the function k(^^^, has the form 

K{x,y) =- + f^^[^{l~ix + iy)-^{l + iy)]]-{y^-y} . (71) 
y i X J 



Firstly, we shall check the accuracy of the approximation of the expression (|67D for cr{uj) 
by the "extended" Drude formula, Eq. (pHI). For that purpose we use the transport spectral 
function a1j.{uj)F{uj) of the form shown in Fig. ||. A comparison between Eqs. (|67D and (|68[) 
is shown in Fig. ^ for constants of coupling Atr = 1, where 

Atr = 2 / — alin)Fin) . (72) 


There are some differences in the low energy regime, u ^ 200 cm^^, shown in Fig. ^, especially 
at low temperatures. One should be careful using the "extended" Drude formula in this 
case. However, even in this region the difference is rather small. The difference between 
the conductivity calculated by the formulae (0) and (^) continues to be small even for a 
constant of coupling A ~ 2. 

The transport relaxation rate 

' ' (73) 



rt*(u;) Ttr{uj) m*{uj) 
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has a universal behaviour for metals with different values of A and different phonon density of 
states. Figure | shows the frequency dependence of the relaxation rate in dimensionless 
form. Here uj^ is the value of the maximum phonon frequency (i. e. the end of the phonon 
spectrum). The six different curves plotted in Fig. |^ are indistinguishable as function of the 
dimensionless energy ujjijj^. The quantity is also rather universal as a function of the 

dimensionless temperature T/vrtUm. Figure |^ shows for T = 10 K and 100 K. Firstly, we 
would like to emphasise the quasi-linear cu-dependence of over a large energy interval, 
0.5 < C(j < (3 — 4) cUin- The function increases with increasing energy uj up to very high 
values, Lo ~ lOc^m. This contrasts strongly with the behaviour of the one-particle relaxation 
rate defined by Eq. (0) and shown in the inset of Fig. |^. It is well known |1^ that the 
latter rapidly increases with energy and becomes constant for a; > lo^. This difference in 
the behaviour of the functions }, ^ and -tW was first discussed in Ref. [Bl. 

As mentioned above, there are only a few investigations of the influence of the 
electron-phonon interaction on the optical spectra of normal state metals where the fre- 
quency dependence of these effects was observed. Besides the reasons mentioned above 
there is another very important reason for the small number of such investigations. Namely, 
the absolute values of the frequency dependence of the discussed effects is determined by 
the value of }, at — 00 [l6|. This function can be written as 

lim — \— = 7i\ (u) , (74) 

where 

00 

X{uj) = 2 Jdio al{uj)F{uj) . (75) 


This value can be expressed as 

^ (1 - 2) A^^ , (76) 

Ttr 

as can be seen from Figs. ^ and ^, and it is very small for ordinary metals. In lead, for 
example, one finds A (u) ^ 100 cm^^. It is extremely difficult to observe this phenomena 
in the usual reflection spectra. Therefore, the observation of Holstein processes in lead was 
made using the light scattering inside a cavity whose walls hold the sample material. 
It was averaged over at least 100 such reflections to increase the accuracy of the experi- 
ment. It was shown |3,|18| that a far-infrared measurement at low temperature, containing 
phonon-induced structure, can be inverted to give the spectral function of the electron- 



phonon interaction a'i^{uj)F(uj). This function was obtained in Ref. |18] for lead in very 
good agreement with experimental data from tunnelling measurements. We shall come back 
to discuss this observation in some more detail. 



III. ELECTRON PHONON INTERACTION AND OPTICAL SPECTRA OF 

HTSC SYSTEMS 

A large amount of work has been devoted to the study of optical spectra of HTSC systems 



see for example the reviews 119,20]). Investigations include normal and superconducting 
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state properties, doping dependence, the effect of impurities, etc. . Here, we shall restrict 
our analysis to optimally doped HTSC in the normal state. In this case the optical spectra 
of all HTSC materials show quite similar behaviour: their reflectivity drops nearly linearly 
with energy from i? ~ 1 down to i? ~ 0.1 at the plasma edge a;*, where the values of u* 



vary for different HTSC materials, 1 eV < cUp < 1.8 eV |2^. Measurements for different 
HTSC compounds further coincide in showing large values for the conductivity in the energy 
interval 2 eV < < 15 eV and well developed charge fluctuation spectra (described by the 
energy loss function — Im (^-^^) in the energy interval 5 eV < < 40 eV The high 
energy part of spectra above a few eV can be well described in terms of the usual band 



structure calculations p0|j2^ . Moreover, the calculations |l22|j23[| can also describe quite 



accurately some low energy interband transitions in YBa2Cu307 

The most unusual part of the optical spectra of HTSC systems is connected with the 
strong decrease in reflectivity, almost linear in energy, in the region < a; < o;* ~ 1 eV [p5 



This behaviour certainly cannot be explained with the simple Drude approach. Thomas 
et al ||2^ proposed two different ways of analysing such spectra. First, a one-component 
Fermi-liquid approach, using the "extended" Drude formula with frequency dependent mass 
m*{ijj) and relaxation rate This model successfully describes the heavy fermion sys- 



tems 12^. Secondly, a two-component approach, where the spectrum is decomposed in two 
components, the Drude part and a mid-infrared (MIR) absorption band. In the latter case 
there is no unique way of separating the two components. 

As we have discussed above, including the electron-phonon interaction in the consid- 
eration leads immediately to the representation of the optical conductivity in terms of the 
"extended" Drude formula. It has been shown earlier for some HTSC that the existence 
of strong electron-phonon coupling, including some amount of MIR excitations, can indeed 
explain their optical conductivity. We now consider this approach in some more detail. To 
describe the optical spectra of HTSC systems we use the formulae obtained in section |I| 
using a transport spectral function al^{uj)F{uj) of the form shown in Fig. |1], multiplied by 
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This spectral function extends up to uj^ = 735 cm~^ and well resembles the phonon spectra 
of HTSC systems. Generally, the optical properties do not depend on the actual shape 
of alj.{u)F{u) but on some moments of it. Thus, we use the constant of electron-phonon 
coupling Atr as deflned in Eq. (|72[) as flt parameter. 

Figure ^ shows the reflectivity of optimally doped La2_xSr2;Cu04 (LSCO) at room tem- 
perature. We have used the values ujp = 1.8 eV and eoo = 4.6 obtained by one of us ||28| from 
band structure calculations. The contribution of the electron-phonon coupling was supposed 
to be Atr = 2.5. The resulting reflectivity corresponds well to measurements on good quality 
films |]23|. The conductivity calculated in Ref. [^] describes well the optical spectra of LSCO 
at energies above ~ 3 eV, as was confirmed by experiment ^9]. Taking all this into account. 



we can conclude that the band structure approach joined with strong electron-phonon in- 
teraction can explain the overall behaviour of the optical spectrum of LSCO in the energy 
range < a; < 40 eV. Moreover, this approach also explains the temperature dependence 
of the optical spectrum of LSCO quite well. Figure |^ shows the refiectivity of LSCO in the 
FIR at temperatures T = 100, 200, 300 K. The agreement between experimental ||30| and 
calculated curves is good considering that the only fit parameter used was the constant of 
coupling Atr. In Fig. we show the refiectivity for LSCO up to u; ~ 1 eV for the same 
temperatures. The agreement between calculated and measured curves is still quite good on 
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this larger energy scale. Nevertheless, we should point out a small discrepancy in the MIR 
region. While this discrepancy is very small at T = 100 K, it becomes more pronounced 
with increasing temperature. Introducing a MIR band, this can be seen as further evidence 
for the temperature dependence of such a band, as discussed in Ref. p. 

We now focus on the temperature dependence of the optical reflectivity in YBa2Cu307_5 
(YBCO), using the values Up = 3 eV and eoo = 5.8 obtained from band structure calculations 
P^ . Figure ^ shows the reflectivity of YBCO from Ref. compared to our calculations 
for T = 100, 200, 300 K, where we have used an impurity scattering rate ^— = 300 

'Hmp 

cm~^ in our calculation of the relaxation rate. The agreement between experimental and 
theoretical curves seems to be again quite good. The only fit parameter used here was 
Atr = 3. Certainly, there are some discrepancies between data and theoretical curves in 
Fig. p. We should point out that discrepancies of nearly the same order exist between 
different type of monocrystals and films. Also, our approach is based on the quasi-isotropic 
approximation for the electron-phonon interaction which is clearly an oversimplification in 
these quasi-two dimensional systems. The inset of Fig. ^ shows the calculated reflectivity 
of YBCO up to 1.8 eV at room temperature. It can be seen that there is an upturn in the 
reflectivity at MIR frequencies which is certainly much larger than in LSCO. This confirms 
the observation made in Ref. [^] that the MIR band is more pronounced in YBCO than 
in LSCO. Furthermore, we know from band structure calculations ||2^ that YBCO systems 
show interband transitions with rather strong intensity in the MIR region. 

Figure |10| shows the behaviour of the transport relaxation rate obtained for YBCO 
at temperatures T = 100, 200, 300 K. It well resembles the usual shape of such curves 



in HTSC and might be compared to Fig. 2 in Ref. [Q, where the relaxation rate for 
Bi2Sr2CaCu208 was derived from the conductivity using the "extended" Drude formula. 
It was also shown there that the behaviour of the relaxation rate ■:pr^ cannot be explained 

in the framework of marginal Fermi liquid theory |^|, since saturates at energies 



uj > ujc — 1500 cm~^ whereas according to marginal Fermi liquid theory the relaxation 
rate should continue in a straight line. This argument was used in Ref. as evidence 
against the one-component approach to the optical spectra of HTSC systems. However, the 
saturating behaviour for 75:^ can be clearly seen in Fig. |TU|, at nearly the same energy as 
it was observed in Bi2Sr2CaCu208. Figure |Tl| shows the calculated resistivity for YBCO, 
demonstrating clearly a linear temperature dependence over a large temperature interval. 

To conclude this considerations, we would like to discuss shortly the possibility of in- 
verting reflectivity data to obtain the transport spectral function a1j.{uj)F{uj). This has 



been done for lead at T = [|18|, where the the infrared data was obtained by applying a 
magnetic field to drive the system into the normal state. There, a spectral function with 
detailed structure could be obtained. In HTSC systems the situation is more complex. The 
superconducting state exists up to considerably high temperatures T > Tc and cannot be 
suppressed by magnetic fields to allow IR measurements at T = 0. At finite temperature 
the inversion procedure is far more difficult and as a result such calculations only yield the 
rough overall form of the transport spectral function. Nevertheless, the resulting a'^j.{uj)F{uj) 
clearly resemble the phonon spectra of these systems M. 
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IV. CONCLUSION 



The results obtained in this work demonstrate clearly, that strong electron-phonon in- 
teraction, combined with band structure calculations, describe the overall behaviour of the 
optical spectra and the main part of the transport properties of HTSC in a straightforward 
manner. However, we do not claim that this simple quasi-isotropic approach can explain 
all details of the behaviour of HTSC systems, even in the normal state, and even at optimal 
doping. There are a number of open problems concerning the behaviour of the Hall coeffi- 
cient, the NMR relaxation rate of the copper sites. We would hke to point out the existence 
of at least two relaxation rates in the HTSC systems, namely the quasi-particle relaxation 
rate and the transport -ipr^ which are very different over a large energy range. We 
cannot rule out that the relaxation rate involved in the Hall current, due to strong and pos- 
sibly anisotropic electron-phonon interaction, will be different from the transport relaxation 
rate in those systems. This could lead to the observed temperature dependence of the Hall 
coefficient. 

Last but not least, we should point out that the simple approach presented here does not 
work at low temperatures. It cannot properly describe the anisotropy of the superconducting 
order parameter, although it yields the correct order of magnitude for the value of Tc. There 
are additional phenomena besides electron-phonon coupling which become important at low 
temperatures. There are a number of different models which combine a strong electron- 
phonon interaction with interband Coulomb interaction, with the existence of a Van Hove 
singularity in the electron spectrum, etc., of which a detailed discussion is beyond the scope 
of this work. 
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FIGURES 



FIG. 1. Transport spectral function a^^{io)F{uj) used in the calculations. 

FIG. 2. Comparison between the calculated optical conductivity (real part) ai{uj) using Eq. 
( pTD (dashed line) and using the "extended" Drude formula, Eq. (^) (solid line). 

FIG. 3. Same as Fig. || at low energies 

FIG. 4. Transport relaxation rate at T = 10 K for different constants of coupling, 

Atr = 0.2, 0.6, 1.0, 1.4, and different phonon spectra with = 612 cm~^, 735 cm~^, 882 
cm~^, calculated using the "extended" Drude formula. In dimensionless units all curves fall on top 
of each other. 

FIG. 5. Transport relaxation rate at T = 10 K and 100 K for A = 1, lo^ = 735 cm~^. 

The inset shows the one-particle relaxation rate calculated for the same parameters. 

FIG. 6. Calculated reflectivity for optimally doped LSCO at T = 300 K, using the "extended" 
Drude formula with A = 2.5, tOp = 1.8 eV, eoo = 4.6, 1/rimp = 100 cm~^. 

FIG. 7. FIR reflectivity for LSCO from Gao et al at different temperatures compared to our 
calculations. 

FIG. 8. Same as Fig. |^ on a larger energy scale. 

FIG. 9. FIR reflectivity for YBCO from Schiitzmann et al at different temperatures compared 
to our calculations, using A = 3, Wp = 3 eV, eoo = 5.8, 1/rimp = 300 cm~^. The inset shows the 
calculated curve at room temperature in a larger energy scale. 

FIG. 10. Calculated transport relaxation rate -prj^ for YBCO at different temperatures. 

FIG. 11. Calculated DC resistivity for YBCO as a function of temperature. 
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Fig. 1 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 
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2 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 




3 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 




Fig. 4 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 




Fig. 5 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 
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6 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 




Fig. 7 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 




Fig. 8 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 
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Fig. 10 Electron-phonon interaction and optical spectra of metals, Kaufmann, Maksimov, Salje 
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